Third-harmonic exponent in three-dimensional iV-vector models 
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Abstract 

We compute the crossover exponent associated with the spin-3 operator in 
three-dimensional O(N) models. A six-loop field-theoretical calculation in the 
fixed-dimension approach gives 03 = 0.601(10) for the experimentally relevant 
case N = 2 (XY model). The corresponding exponent (3% = 1.413(10) is 
compared with the experimental estimates obtained in materials undergoing 
a normal-incommensurate structural transition and in liquid crystals at the 
smectic-A-hexatic-B phase transition, finding good agreement. 

PACS Numbers: 05.70.Jk, 64.60.Fr, 75.40. Cx, 61.30.-v 
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In nature many physical systems undergo phase transitions belonging to the universality 
classes of the 0(N) vector models. In particular, the XY model, corresponding to N = 2, 
describes the A-transition in 4 He, (anti)ferromagnets with easy-plane anisotropy, density- 
wave systems, etc.; see Ref. [1] for a review. The critical exponents associated with the 
order parameter have been accurately measured both experimentally and theoretically [1]. 
Moreover, in some XY systems it is also possible to measure experimentally the critical 
exponents associated with secondary order parameters. This is the case of liquid crystals 
[2-4], of normal-incommensurate transitions [5-7], and of graphite- intercalation compounds 
[8]. 

The most relevant exponent is the second-harmonic one that has been recently computed 
to high-precision using field-theoretical methods in Refs. [9,10]. The e-expansion gives 02 = 
1.174(12), while the fixed-dimension expansion gives 02 = 1.184(12). The fourth-harmonic 
crossover exponent was reported in Refs. [11,10]: 4 = —0.077(3) (e expansion) and 4 = 
—0.069(5) (fixed-dimension expansion). Here, we wish to determine the third-harmonic 
exponent by means of a six-loop perturbative calculation in the fixed-dimension approach, 
extending previous three-loop determinations [12,13,3]. Such a calculation is also relevant 
for some crossover phenomena, in which the XY symmetry is reduced to that of the three- 
state Potts model, as it happens in cubic magnets in the presence of stress or of appropriate 
magnetic fields [14,15]. 

In the field-theoretical approach one starts from the usual 4 Hamiltonian 



n = J d d x [^0) 2 + ir0 2 + ^(0 2 ^ 2 



(1) 



where <j> a (x) is an iV-component real field. The XY model corresponds to N = 2, but here 
we will keep N generic. Secondary order parameters are associated with operators 
that are polynomials of order / in the fields and that transform irreducibly under the spin-/ 
representation of the O(N) group. In particular, the third-harmonic operator is 

3 2 
Oabc = <Pa<Pb<j)c ~ N + 2 (0Ac + 4>b$ac + 4>c$ab) • (2) 

We wish now to compute the crossover exponent 3 associated to 0^ c and the corresponding 
exponents fi 3 and 73 given by 

P 3 = 2 - a - 03, 

73 = -2 + a + 20 3 . (3) 

The exponents (3$ and 73 describe respectively the critical (singular) behavior of the average 
(O^(x)) ~ and of the susceptibility xo = E x (O {3) (0)O^(x)) c ~ \t\~^. 

For this purpose we determine the renormalization function Z 3 (g) from the one-particle 
irreducible three-point function (0) with an insertion of the operator 0^ c at zero external 
momenta, i.e. we set 

(0^ a M c ) lpI = AZ^(g), (4) 

where A is a numerical coefficient that ensures that Z 3 (0) = 1, and g is the four-point 
renormalized coupling. Then, we compute the renormalization-group function 
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and 773 = 773(5*), where g* is the fixed-point value of 5. Finally, the renormalizat ion-group 
scaling relation (valid in three dimensions) 

03 = (rj3 + I - ^) f (6) 

allows us to determine 3 . 

We computed T 3 ; (0) to six loops. The calculation is rather cumbersome since it requires 
the evaluation of a few thousand Feynman diagrams. We handled it with a symbolic ma- 
nipulation program, which generates the diagrams and computes the symmetry and group 
factors of each of them. We used the numerical results compiled in Ref. [16] for the integrals 
associated with each diagram. We obtained 

6 _ 2(iV + 10)_ 2 128.736 + 15.4900 iV — 0.650238 N 2 _ 3 
- ( iV + 8 ) 2 9 (iV + 8) 3 9 

1148.68 + 191.005 N + 1.82163 N 2 + 0.283028 N 3 _ 4 
+ (A^ + 8) 4 9 

12606.9 + 2550.46 N + 64.4818 N 2 - 2.34060 N 3 - 0. 152501 N 4 _ 5 

(iV + 8) 5 9 

161373. + 38736.8 N + 1874.23 N 2 - 5.98451 N 3 + 1.88168 A^ 4 + 0.094179 N 5 _ 6 
+ (N + 8) 6 9 

+0(g 7 ), (7) 

where, as usual, we have introduced the rescaled coupling g defined by 

Here g is the usual four-point renormalized coupling normalized so that g = u/m (m is the 
renormalized mass) at tree level. 

Field-theoretical perturbative expansions are divergent, and thus, in order to obtain 
accurate results, an appropriate resummation is required. We use the method of Ref. [17] 
that takes into account the large-order behavior of the perturbative expansion, see, e.g., 
Ref. [18]. Mean values and error bars are computed using the algorithm of Ref. [11]. 

Given the expansion of 773(g), we determine the perturbative expansion of 03(g), 03(g), 
and 73(g) using relations (3) and (6). Then, we resum the perturbative series and compute 
them at g = g* [19]. For N = 2 we obtain 3 = 0.5963(21), 0.5968(2), = 1.398(8), 
1.405(3), 73 = —0.800(7), —0.808(13), where for each exponent we report the estimate 
obtained from the direct analysis and from the analysis of the series of the inverse. The 
two estimates obtained for each exponent agree within error bars, but, with the quoted 
errors, the scaling relations (3) are not well satisfied. For instance, using v = 0.67155(27) 
(Ref. [20]) and /3 3 = 1.403(8) we obtain 3 = 0.611(8), while using the same value of v 
and 73 = —0.803(13) we have 3 = 0.606(6). These two estimates are slightly higher than 
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TABLE I. Critical exponents associated with the spin-3 operator O 
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those obtained from the analysis of 4>z{g) and l/<f> 3 (g). Clearly, the errors are somewhat 
underestimated, a phenomenon that is probably connected with the nonanalyticity [21-23] 
of the renormalization-group functions at the fixed point g*. 

In order to obtain a conservative estimate, we have thus decided to take as estimate of 
03 the weighted average of the direct estimates and of the estimates obtained using fa and 
73 together with the scaling relations (3). The error is such to include all estimates. The 
other exponents are dealt with analogously [24]. The final results for several values of N are 
reported in Table I. 

Let us compare these results with previous ones for N = 2. Ref. [3] reports fa = 
3/? + 3z/x 3 , where ~ 0.174 from the e-expansion and x 3 0.276 from the fixed- dimension 
expansion. It follows fa ~ 1.397 and fa m 1.602 in the two cases. These results are rea- 
sonably close to ours. The exponent fa has been determined at the smectic-A-hexatic-B 
phase transition in liquid crystals, obtaining [2,3] fa ~ 4.8(3 ~ 1.66, using (3 = 0.3485(4) 
(Ref. [20]). Such an exponent has also been measured in some materials exhibiting a struc- 
tural normal-incommensurate phase transition. From the analysis of x-ray scattering data 
in erbium Ref. [6] obtains fa = 1.8(3), while two different experiments in Rb2ZnCl4 give 
respectively fa = 1.80(5) (Ref. [5]) and fa = 1.50(4) (Ref. [7]). Keeping into account that 
the experimental errors seems to be underestimated, there is reasonable agreement with our 
results. The exponent 5* = (3/(f) 3 was measured at the trigonal-to-pseudotetragonal tran- 
sition in [lll]-stressed SrTi0 3 obtaining [14] 5* = 0.62(10). Such an estimate is in good 
agreement with our prediction 5* = 0.58(1). Finally, we should mention that our results are 
also relevant for polymer physics. Indeed, we can derive from the estimates obtained for 
N = the partition-function exponent p for nonuniform [25] star polymers with three arms: 
p = 3(7 + v)/2 + 03 = 3.06(1), where we used 7 = 1.1575(6) (Ref. [26]). 
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